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CHARTS RELATING THE COMPRESSIVE BUCKLING STRESS OF LONGITUDINALLY 
SUPPORTED PLATES TO THE EFFECTIVE DEFLECTIONAL AND 
ROTATIONAL STIFFNESS OF THE SUPPORTS ^ 


B3' Rocier a . Anderson and Joseph W. Semonian 


SUMMARY 

A stability analysis is made of a long flat rectangular plate 
subjected to a uniform longitudinal compressive stress awl 
supported along its longitudinal edges and along one or more 
longitudinal lines by elastic line supports. The elastic supports 
possess deflectional and rotational stiffness. Such a configura- 
tion is an idealization of the compression cover skin and internal 
structure of wing and tail surfaces. The results of the analysis 
are presented in the form of charts in lohich the buckling-stress 
coefficient is plotted against the buckle length of the plate for 
a wide range of support stiffnesses. The charts make possible 
the determination of the compressive buckling stress of plates 
supported by members whose stiffness may or may not he 
defined by elementary beam bending and twisting theory but 
yet whose effective restraint is amenable to evaluation. The 
deflectional and rotational stiffness provided by longitudinal 
stiffeners and full-depth webs is discussed and numerical 
examples are given to illustrate the application of the charts 
to the design of wing structures. 

INTRODUCTION 

In current thin-wing construction, thick cover skins are 
often supported or stiffened by thinner gage internal members 
whose stiffness determines the stabilitj^ and strength of the 
cover skins. A careful evaluation of this stiffness is required 
for members such as longitudinal stringers and full-depth 
webs whose behavior may be substantially influenced by 
local bending of riveted attachment flanges and by shearing 
deflections. When such distortions are present, cover-skin 
budding stresses are usuallj^ overestimated by the usual 
stability criteria which are based upon idealizations of the 
supporting members as beams (or plates) integrally joined 
to the cover skin and possessing stiffnesses El and GJ defined 
by elementary bending and twisting theory. This is borne 
out by a number of tests — for example, references 1 to 3 — 
in which large reductions in budding stress (and failing 
stress) from theoretical values based dn integral support 
theories are reported. The desirability of relating plate 
stability to a stiffness parameter which defines the actual 
or effective stiffness provided by supporting members is 
therefore evident. 


Reference 4 describes a mode of instability of cover skins, 
denoted as wrinlding, the occurrence of which in skin 
stringer panels is attributed to flexibility of the attachment 
flanges of the stringers. In reference 5, this same mode is 
described and is called forced crippling. An approximate 
stability analysis which takes into account flange flexibility 
is given for plates supported by longitudinal stringers or 
by full-depth webs as in a multiweb wing. 

The purpose of this report is to present stability criteria 
which apply to the wrinkling as well as to the more usual 
local instability modes for a number of supported plate 
configm’ations frequently occurring in aircraft-wing con- 
struction. In the design charts presented, the elastic- 
buckling-stress coefficient is given as a function of the buckle 
length of the cover skin for the practical range of effective 
deflectional or torsional stiffnesses of supporting members. 
A section of the report is devoted to a discussion of procedures 
for evaluating the effective deflectional and torsional stiffness 
provided by longitudinal stringers and full-depth webs. 
Numerical examples are then given which illustrate this 
evaluation for practical design cases. The derivations of 
the stability criteria are included in the appendixes. 

SYMBOLS 

b width of plate between intermediate supports 

X length of buckles 

(3=X/6 

t thicloiess of plate 

x,y coordinate axes in length and width direc- 

tions, respectively 

w deflection normal to plane of plate 

p number of bays in width of plate 

2 number of bucldes occurring across width of 

plate 

n integer 

Fourier coefficients 

N compressive load per unit width acting in 

x-direction (length direction) required to 
cause buckling 

k nondimensional buckling-load coefficient, 

Nbf irW 

a compressive stress 

(Ter critical compressive stress 
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Young’s modulus of elasticity 
Poisson’s ratio 

plate flexural stiffness per unit width, 

Et^ 

12 ( 1 -/) 

deflectional stiffness per unit length of support, 
Ib/in.^ 

rotational stiffness of intermediate support 
(moment per unit length required to pro- 
duce a rotation of 1 radian) 
rotational stiffness of edge support (moment 
per unit length required to produce a 
rotation of 1 radian) 

nondimensional deflectional restiaint param- 
eter 

nondimensional rotational restraint param- 
eters 

nondimensional rotational restraint parameter 
fi’om reference 6 

plate edge rotational stiffnesses defined in 
reference 7 

plate carryover factor defined in reference 7 
energies of deformation 
work of applied stress 
total potential energy of system 
energjr parameter 
Lagrangian multipliers 

coefRcients defining amplitude of support 
deflection 

cross-sectional area of stiffener 
moment of inertia of stiffener cross section 
about its own center of gravity 
modal coefficient affecting deflectional stifl'- 
ness of longitudinal stiffener 
nondimensional bending stifl’ness parameter 
for stiffeners of sturdy cross section 
ratio of average stress in stiffener to average 
stress in plate 
Euler column load 
torsion constant 
shear modulus of elasticity 
torsion coefficient which takes into account 
bending stiffness 
polar moment of inertia 

amplitude of sinusoidally distributed lateral 
load 

lateral deflection of longitudinally compressed 
stiffener subjected to sinusoidal lateral load 
depth of web 
thickness of web 

plate flexural stiffness per unit width of web, 
EU/ 


kw budding-stress coeffident of web 

z distance between center of gravity of stiffener 

and middle plane of plate 

p radius of gyration of stiffener about its 

centroid 

SCOPE OF ANALYSIS 

In figure 1 are shown portions of several wing cross sections 
in which the material carrying bending stress is mainly con- 
centrated in the thick plates forming the wing contour. 
Running spanwise are a number of lighter structural mem- 
bers in the form of longitudinal stiffeners and full-depth 
webs. In addition to carrying longitudinal stresses these 
members resist cover-plate deflection and rotation at their 
respective locations by virtue of their stiffness. If the stiff- 
ness characteristics of these members can be defined, the 
buckling stress for the construction can be calculated. 

Ill this analysis the assumption is made that longitudinal 
stiffeners and full-depth webs will pi’ovide a restraint to the 
attached cover plate which is proportional to the distortions 
of these support members. This condition is met if sinus- 
oidally distributed normal loads or torsional moments on 
the supports are assumed to cause sinusoidally distributed 
distortions which are in phase with the loading. Thus sup- 
port stiffness, which is the ratio of load intensitj^ to distortion 
at any point, is a constant along the length of the support. 
With this support characteristic, the attached plate will 
biiclde with deflections and rotations that are distributed 
sinusoidally in the length direction. 

A cross section of the cover-plate buckling modes consid- 
ered most likely to occur are sketched at the right of each 
wing cross section in figure 1 and are denoted cases 1 to 6. 
Cases 1, 2, and 3 primarily involve the deflectional stiffness 
chai’actci'istics of the support members, and cases 4, 5, and 6 
involve the torsional stiffness characteristics of the supports. 
For a given wing cross section, both modes of budding 
shoidd be investigated to determine which mode leads to the 
lower budding stress. 

Cases 1 and 4 represent the buckling modes of a cover 
plate supported by substantial shear webs with an inter- 
mediate spanwise member (shown as a longitudinal stiffener) 
centrally located between the webs. The shear webs are 
assumed to prevent deflection but ma}'- offer a torsional 
restraint to the cover plate. In case 1 the stabilit}’^ of the 
compressed plate was investigated for a range of deflectional 
stiffnesses of the intermediate support and in case 4 the 
torsional stiffness of the supports was considered. Because 
the two lowest buckling modes arc either symmetrical or 
antisjmmetrical with respect to the spanwise center line of 
the plate, it is not necessary to consider both the deflectional 
and rotational stiffnesses of the support simultaneously. 

Cases 2 and 5 represent the most likely budding modes 
for a cover plate rvith two equally spaced spanwise stiffening 
members of equal stiffness between shear webs. In case 2 
the effect of support deflectional stiffness was investigated 
by assuming the torsional stiffness of the intermediate sup- 
ports to be zero. The torsional stiffness of the intermediate 
supports was considered in case 5 with the assumption that 
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Cross sections 



Buckling modes considered 






FirjuRE 1. — Thick-skin box-boani cross sections and the buckling modes considered for each. 


the supports are capable of preventing plate deflection at 
their locations. 

Cases 3 and 6 represent the most likely buckling modes 
for a plate stabilized by many spanwise lines of support of 
identical stiffness. These supports may be full-depth webs, 
as indicated in figure 1, or longitudinal stiffeners. In case 3, 
the deflectional stiffness of the supports was considered by 
assuming the support torsional stiffness to be zero. Toi- 
sional stiffness of the supports was considered in case 6 in 
which the deflections along the supports arc assumed to be 
zero. 

The loading and support conditions for the si.x cases con- 
sidered are shown schematically in figure 2. 'Fhe com- 
pression cover plate is represented by a uniformly compressed 
long flat plate which is simply suppoi'tcd at the loaded edges. 
The deflectional stiffness of the supports is represented by 
an elastic spring whose stiffness per unit length is denoted 
by The stiffness ip may include the flexibility of the 
tension cover of a midtipost stiffened wing (ref. 8) in which 
tension cover flexibility would have an effect on the stability 
of the compression cover. The parameter ip as defined in 
this report is a generalization of the foundation modulus 
concept as used by Timoshenko for beams on an elastic 
foundation (ref. 9). The support torsional stiffness param- 
eters are denoted by y and a. The parameter y is associated 
with the torsional stiffness of the nondeflecting shear webs 
and a is associated with the torsional stiffness of the inter- 
mediate supports. These two parameters are equivalent to 
the torsional stiffness parameter 4So defined by .Lundquist 
and Stowell in reference 6. 

For each of the first three cases a stability criterion in 
closed form is derived by the Lagrangian multiplier method 
(ref. 10). For the last three cases a stability criterion is 


ol)taincd I),v using the principles of moment distribution 
explained in reference 7. With these stability criteria, 
numerical calculations have been made and are presented 
in design-chart form. 

PRESENTATION OF STABILITY CRITERIA 

Cases 1, 2, and 3. — The stability criteria for cases 1, 2, 
find 3 which involve the deflectional stiffness of the inter- 
mediate supports are presented in appendix A as equations 
(Al9), (A24), and (A28). In these equations, the effective 
deflectional stiffness \p of the supports is contained in the 
nondimensional parameter \pPlir*D, and the effective tor- 
sional stiffness y provided along the shear webs is contained 
in the nondimensional parameter yblir^D. Values of the 
parameter ipPlir^D may be determined from these equations 
as a function of the compressive-buckling-stress coefficient 

and the ratio of buckle length to bav width \/b for 

w U 

assigned values of the torsional restraint parameter yhlw^D. 

Two sets of numerical calculations have been made by 
assigning the values 0 and to yhjiv^D-, these values corre- 
spond to simple support and complete fixity, respectively, 
along the shear webs. These numerical results are presented 
in tables I, II, and III. Cross plots of the values in the 
tables have been made to form design charts (figs. 3 to 7). 
From these charts^, the combinations of ipb^jir'^D, k, and \/b at 
which buckling is initiated, may be read. The cutoffs in 
figures 5 and 7 define the values of iPP/tt^D at which general 
instability involving deflection of the cover and the supports 
changes to local buckling of the cover (no support deflection) 
in accordance with the assumption made that the supports 
possess zero torsional stiffness. 
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Figure 2. — Six cases for which stability criteria are presented. 


In oi'dor to use the charts for plates on particular types of 
supports, the parameter \p¥lir*D for the support must be 
evaluated. For the usual type of support, such as a longi- 
tudinal stiffener, or a full-depth web, will be a 

function of the stresses in the support and the wave length 
of buckling, as well as the physical characteristics of the 
support. A discussion of the evaluation of for longi- 

tudinal stiffeners and webs is given in the section entitled 
“Effective Stiffness of Supports,” and numerical examples 
illustrating the procedure are given in a subsequent section 
entitled “Illustrative Examples.” 

Cases 4 , 5 , and 6. — ^For cases 4 , 5 , and 6, the cover is 
restrained by equally spaced nondcflecting supports of equal 
rotation stiffness a while the plate side edges are restrained 
by nondeflecting supports of equal rotational stiffness 7. 
The stability criteria for these cases are given in appendix B 
as equations (.B2), (B6), and (BlO). Values of the rotational 
stiffness parameter ablwW required to develop a given com- 


pressive-buckling-s tress coefficient k= 


Nb^ 

ttW 


in the cover at a 


given ratio of buckle length to bay width X/6 may bo deter- 
mined from these equations for assigned values of the edge- 
restraint parameter yblw^D. As was done for the deflectional 
stiffness cases, numerical results are presented for 76/ir^Zt 
equal to 0 a.nd 00 . The numerical results were obtained by 
using the stiffness tables of reference 11 and have been 
plotted to form design charts (figs. 8 to 12). 

For a given design problem in which the supports have both 
deflectional and rotational stiffness, the buckling-stress coeffi- 
cient obtained by considering the mode of buckling which 
involves the rotational stiffness of the supports must be 
compared with the coefficient obtained by considering the 
mode involving primarily the deflectional stiffness of the 
supports. The lower of these two values defines the buckling 
stress for the configuration. The evaluation of the torsional 
stiffness of longitudinal stiffeners and full-depth webs is 
discussed in the next section. 
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Figurk 3. — Stability curves for case 1 with simply supported side edges. 


J ^ 35 / 2.5 /, 8 /|. 4 /|I .0 /.7/.5 / .3 / .1 ,/ o / -.08 /-.I2 


MIlillB 


8 9 







COMPRESSIVE BUCKLING OF ELASTICALLY SUPPORTED PLATES 


7 



Figure 7. — Stability curves for case 3. 


TABLE I 


VALUES OF DEFLECTIONAL-STIFFNESS PARAMETER 

FOR CASE 1 





yph^jir^D for values of k 



X/6 








0 


2 

3 

4 

5 

6 


"yf) 

Simply supported side edges; 

0.4 

-19. 89 

-17. 85 

-15. 67 

-13. 32 

-10.71 

-7. 690 

-3. 928 

.5 

-10. 19 

-8. 526 

-6. 680 

-4. 534 

-1.827 

2. 253 

12. 26 

.6 

-5. 899 

-4. 497 

-2. 861 

-.7810 

2. 376 

9. 977 


.7 

-3. 721 

-2.517 

-1.050 

.9544 

4. 658 

18.79 


.8 

-2. 504 

-1.401 

-. 1573 

1.709 

5. 309 

23. 10 


.9 

-1. 773 

-.8683 

.2748 

1.927 

5. 076 

18. 52 


1.0 

-1.309 

-. 5230 

.4676 

1.877 

4. 411 

12. 58 


1.2 

-.7887 

-. 1916 

.5426 

1.527 

3. 052 

6.177 

20. 15 

1.6 

-. 3826 

-.0197 

.3999 

.9039 

1.542 

2.411 

3.760 

1.8 

-.2948 

-.0034 

.3254 

.7062 

1.162 

1.733 

2.495 

2.0 

-. 2384 

0 

.2638 

5610 

.9032 

1.309 

1.807 

2.5 

-. 1624 

-.0081 

.1571 

.3353 

.5293 

.7427 

.9803 

3.0 

1268 

-.0192 

.0937 

.2127 

.3387 

.4725 

.6157 

4.0 

-.0956 

-.0349 

.0275 

.0917 

.1580 

.2263 

.2970 

6.0 

-.0754 

-.0488 

-.0215 

.0062 

.0342 

.0626 

.0914 

8.0 

-.0688 

-.0542 

-.0389 

-.0235 

-.0080 

.0076 

.0233 

10.0 

-.0672 

-.0568 

-.0471 

-.0373 

-.0274 

-.0175 

-.0076 

Clamped side edges: « 

■k^D 

0.4 

-19. 895 

-17.85 

-16. 67 

-13.33 

-10. 73 

-7. 766 

-4.155 

.5 


-8.541 

-6. 717 

-4.626 

-2. 092 

1.338 

7.154 

.6 

-5.9169 

-4. 5417 

-2.960 

-1.079 

1.417 

5. 371 

14.86 

.7 


-2. 606 

-1.266 

.3898 

2. 658 

6. 424 

16.11 

.8 

-2. 565 

-1.596 

-.4663 

.9286 

2. 818 

5. 828 

12. 499 

.9 


-1.041 

-.0965 

1.0513 

2. 555 

4. 777 

8.863 

1.0 

-1.4087 

-, 7233 

.0654 

1.004 

2. 186 

3. 810 

6. 376 

1.2 


-.4208 

. 1383 

. 7765 

1.529 

2. 455 

3. 667 

1.4 


-.3039 

. 1060 

.5594 

1.070 

1.659 

2. 358 

1.6 


-.2655 

.0559 

.3924 

.7602 

1. 168 

1.626 

2.0 

-.4133 

-.2262 

-. 0304 

. 1754 

.3928 

.6234 

.8696 

2.5 


-. 2231 

-.0997 

.0276 

. 1593 

.2957 

.4374 

3.0 


-. 2263 

-.1416 

-.0549 

.0338 

. 1246 

.2177 

4.0 

-.2796 

-.2330 

1859 

-. 1381 

-.0897 

-.0408 

.0088 

6.0 

-.2604 

-.2398 

-.2190 

-. 1982 

-. 1771 

-.1560 

1347 

8.0 

—.2542 

—.2426 

-.2308 

2186 

-.2074 

-. 1958 

-.1839 

10.0 

-.2506 

-.2438 

-.2365 

-.2289 

-.2215 

-. 2159 

-.2065 


TABLE II 


VALUES OF DEFLECTIONAL-STIFFNESS PARAMETER 
FOR CASE 2 





■^b'^lTT^D for values of k 



X//) 








0 


2 

3 

4 

5 

6 


Simply supported side edges; ■;^=0 

0.4 

-19. 826 

-17. 744 

-15. 491 

-13. 029 

-10. 234 

-6. 927 

-2. 654 

.5 

-10.050 







.6 

-5. 707 

-4. 1853 

-2. 368 

-.0123 

3. 511 

10. 600 


.8 

-2. 275 

-1. 1166 

.3343 

2. 345 

5. 626 



1.0 

-1.094 

2287 

.8404 

2. 263 

4. 4295 



1.2 


. 04579 

.8196 

1. 790 

3. 094 

19. 46 


1.4 


. 1263 

.6970 

1.377 

2. 223 

4.82 


1.6 


. 1417 

.5747 

1.0716 

1.656 



2.0 

1279 

. 1218 

.3924 

.6882 

1.016 

1.381 

1.796 

2.5 


.0877 

.2567 

.4366 

.6261 

.8297 

1.048 

3.0 

-.04899 

.0624 

. 1780 

.2982 

.4236 

. 5545 

.6915 

4.0 


. 0325 

.0964 

.1597 

.2284 

.2968 

.3667 

6.0 


.0085 

.0365 

.0648 

.0933 

. 1221 

.1513 

8.0 


-. 000436 

.0162 

.0310 

.0468 

.0628 

.0780 

10.0 

-.01463 

-.00466 

. 00534 

.0154 

.0254 

. 0356 

.0457 



Clamped side edges: 

■K^JJ 

= » 



0.4 

-19.023 

-17.745 

-15. 493 

-13. 035 

-10. 248 

-6. 963 

-2. 758 

.6 


-8. 348 

-6. 327 

-3. 980 

—.9679 

3. 466 

12.90 

.6 

-5. 5807 

-4. 206 

-2.418 

-.1383 

3.178 

10.118 


.8 

-2. 1971 

-1.173 

.2171 

2. 107 

6. 371 



I.O 

-1.0593 

-.3041 

.7094 

2. 0576 




1.2 

-.5948 

-.0334 

. 70025 

1.630 

3. 031 



1.4 

~. 3772 

.0494 

.5919 

1.247 

2.122 



1.6 

-.2637 

.0686 

.4813 

.9603 

1. 550 



2.0 

-. 1602 

.0554 

.3143 

. 5995 

.9209 

1. 296 

1.767 

2.5 


.0268 

. 1890 

.3614 

.5464 

. 7475 

.»9689 

3.0 


.00478 

. 1159 

.2317 

.3530 

.4807 

.6158 

4.0 

-.07566 

-.0215 

.0399 

.1030 

.1673 

.2334 

.3013 

0.0 

-.06696 

-.0430 

-.0160 

.0113 

.0388 

.0666 • 

.0947 

8.0 


-. 0513 

-. 0358 

-.0207 

-.00540 

.00996 

.0243 

10.0 

06336 

-.0544 

-.0451 

-.0353 

-.0257 

-.0159 

-.00619 
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TABLE III 


VALUES OF DEFLECTIONAL-STIFFNESS PARAMETER 

FOR CASE 3 


X/6 

^f/byir*D for values of k 


U 


2 

3 

4 

5 

6 

0.4 

-10. 758 

-17. 622 

-15. 304 

-12. 730 

-9. 775 

-6. 197 

-1. 479 

.5 

-9. 915 

-8. 076 

-5. 955 

-3. 381 

0 

5. 072 

14. 795 

.6 

-5. 524 

-3. 895 

-1.922 

.6415 

4. 366 

10. 933 


.7 

-3. 297 

-1.861 

-.0828 

2. 291 

5.851 



.8 

-2. 064 

-. 8232 

.7235 

2. 720 

5. 915 



1.0 

-.9180 

0 

1. 1066 

2. 502 




1.2 

-.4600 

.2169 

1.004 

1.957 

3. 104 



1.4 


.2570 

.8317 

1.490 

2. 260 



2.0 

-.06203 

. 1919 

.4626 

. 7525 

1. 065 

1. 402 


2.5 


. 1367 

.3059 

.4827 

.6678 

.8621 

1.0664 

3.0 

-.01234 

. 1000 

.2158 

.3351 

.4582 

.5854 

.7169 

4.0 

-. 003909 

.05905 

. 1231 

. 1882 

.2544 

.3219 

.3905 

5.0 


. 03860 

. 07923 

. 1203 

. 1618 

.2038 

.2463 

6.0 

-.000769 

. 02710 

.05519 

. 08349 

.1120 

. 1407 

. 1697 

8.0 

-.000243 

.01541 

.03114 

. 04693 

. 06278 

. 07871 

. 09470 

10.0 

-.00010 

. 009913 

.01995 

.03002 

.04012 

.05024 

. 06039 



\/d 


PiGUHE 8. — Stability curves for case 4 with simply supported side 

edges. 



EFFECTIVE STIFFNESS OF SUPPORTS 

General design charts have been presented, which, with 
one reservation, are independent of the medium providing 
restraint to the compression plate. The reservation is that 
the supporting medium must be of such a type that sinusoid- 
ally distributed normal loads and torsional moments cause 
sinusoidally distributed distortions which are in phase with 
the loading. Such behavior is characteristic of beam stiff- 
ness, as provided by longitudinal stiffeners of sturdy cross 
section. The buckling distortions of the webs of a multiweb 
beam also appear to be distributed sinusoidally along the 
length of the beam, and the reactions of the attachment 
flange on the compression cover of the beam are assumed to 
be proportional to the distortions. 

The inclusion of the effects of cross-sectional distortion 
and shear distortion in the evaluation of the stiffness para- 
meters \ph^lir*D and ablir^D for these two types of supports is 
discussed below. 
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\/b 

Figure 10. — Stability curves for case 5 with simply supported side 

edges. 


Stiffness of longitudinal stiffeners. — The most common 
type of supporting medium for plates is the longitudinal 
stiffener which participates in carrying the compressive load. 
If the distortion characteristics of such a stiffener are defined 
by elementary beam bending theory, the deflection under a 
lateral load of amplitude ;j distributed sinusoidally over a 
length X is 


(A* 


S(x)- 


tt^EL 


sin 


«// 


(xA-s 







deflection, then is 


If the average stress a in the stiffener is proportional to the 
compressive buclding stress acting in the attached plate, 
^ may be written as 


or • 







ck 


bt 6^ 


( 2 ) 


where cAs is the end load carried by the stiffener, and le// is 
the moment of inertia of the stiffener cross section about an 
axis lying in a plane parallel to the attached plate. The 
stiffness of the stiffener, defined as the ratio of lateral load to 


where c is the ratio of the average stress in the stiffener to 
the average stress in the plate. 

The theoretical analysis of reference 12 shows that the effec- 
tive moment of inertia of longitudinal stiffeners attached to 
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one side of a uniformly compressed plate may be expressed 
as a correction to the moment of inertia of the stiffener about 
its own center of gravity Is- In this form, equation (4) of 
reference 12 may be written 


EIs,f_EIs 

bD bD 


1 



(•■ 1 ) 


In equation (3) the quantity 2 is the distance between the 
center of gravity of the stiffener and the middle plane of the 
plate, and p is the radius of g 3 Tation of the stiffener. The 
modal coefficient is a function of budding mode and 
associated wavelength. The variation of with \/b taken 
from reference 12 (which is applicable when the plate side 
edges are simply supported and when Poisson’s ratio is %) 
is given in figure 13. The subscript p denotes the number of 
bays in the width of the plate, and q denotes the number of 
buckles across the width of the plate {q is equal to 1 for the 
cases considered in this report). With EEfflbD defined, 
equation ( 2 ) should give satisfactorv values of the stiffness 
parameter yl/¥/ir*D for stiffeners of sturdy' cross section; 
that is, stiffeners whose cross-sectional and shearing distor- 



Fkjure 13. — Functions appearing in expression for effective flexural 
stiffness of stiffeners attached to one side of plate (from ref. 12). 

(a) Loads on web. 

(b) Idealization of web. 

(c) Deformed shape of idealized web. 


tions under load introduce deflections which are small com- 
pared with the overall deflection as a beam. 

In practical applications stiffeners are often formed from 
sheet, which necessitates a bend radius between the web of 
the stiffener and the attachment flange. For certain propor- 
tions, deflection of the plate may be appreciably increased by 
the flexibility of the attachment flange between the rivet line 
and the web of the stiffener and by shearing distortion in the 
stiffener. If the total deflection 6 is assumed equal to 
61 + 52+ §3 where 61 is the deflection due to bending of the 
stiffener as a beam, 82 is the deflection due to flexibility of the 
stiffener attachment flange, and 63 is the deflection due to 
shearing distortion in the stiffener, the effective stiffness may 
be written as 


-=—+—+— 

^ 'Pl p2 p3 

In nondimensional form the effective stiffness is given b>' 


■ipb^ 1 

tt'^D Tr'i) . tt'^D . tt^D 
pib^ 4/2b^ Psb^ 


(4) 


where piPlir*D is given by the right-hand side of equation ( 2 ), 
must be evaluated either analytically or experimen- 
tally, and may be calculated. It is evident that if 

either pi, p 2 , or 1^3 approaches zero, the effective stiffness of 
the stiffener approaches zero. Any other significant 
distortions can be included in a similar manner. . 

The torsional restraint furnished a plate by a stiffener 
which undergoes no cross-sectional distortion when it twists 
is discussed in reference 13. The expression for its stiffness 
(eq. ( 8 ) of ref. 13 rewritten in the notation of the present 
report) is 
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where the quantities J, Cbt, and Ip must be calculated with 
respect to an assumed axis of rotation. In nondimensional 
form, the stiffness is 



¥ ECbt 


ckir^ 



(5) 


Expressions similar to equation (5) should be derived for 
those stiffeners in which torsional moments applied to the 
stiffener attachment flange cause distortion of the cross 
section of the stiffener when it twists. 

Stiffness of full-depth webs. — When the compression cover 
of a beam is supported by full-depth webs as in a multiweb 
beam, the effective stiffness of the webs in resisting sheet 
deflection and rotation at the skin-web juncture must be 
evaluated. Keference 14, for example, evaluates the effec- 
tiveness of integrally joined webs as torsional restraints on 
the cover of a multiweb beam. The assumption made in 
that analysis is that the webs possess sufficient deflectional 
stiffness to form longitudinal nodes along the skin-web 
juncture during buckling. The range of skin and web 
proportions for which this assumption is valid, however, is 
not established. 

For built-up construction, the deflectional stiffness pro- 
vided by an unstiffened web plate is influenced by the eccen- 
tricity of the connection between web and cover plates and 
by the state of stress existing in the webs of a beam under- 
load. In particular, for channel-type webs formed from 
sheet, appreciable distortions of the attachment flanges and 
lateral deflection of the web are produced by either depthwise 
crushing or stretching forces. In accordance with the 
stifl'ness analysis for longitudinal stiffeners, the stiffness of the 
channel should be analyzed under the action of a depthwise 
load applied sinusoidally along the length of the attachment 
flange in the presence of the stresses that exist in the web 
during beam bending. This procedure is illustrated by a 
numerical example in the next section. The outcome of such 
an analysis is influenced rather strongly by the assumed 
eccentricity of the applied load (with respect to the plane 
of the web) and by the degree of clamping that is assumed to 
be provided by the riveted connection between web attach- 
ment flanges and the cover plates. The importance of these 
factors in calculating deflectional stiffnesses has been 
emphasized in reference 5. 

With regard to the torsional restraint provided to the com- 
pression cover by integrally joined webs, the restraint data 
presented in figure 9 of reference 14 are analagous to equation 
(5) for the torsional stiffness of a stiffener; that is, the 
restraint coefficient e in figure 9 of reference 14 is a measure 
of the negative of the stiffness of a web subjected to a pure 
bending stress distribution as a function of buckle length. 
The relationship between the torsional stiffness parameter 
ablir^D of the present report and the restraint coefficient e is 


ab _ e b Dw 

^bw D 

Wlien webs are not integrally joined to the cover, the stiffness 
of the attachment should be taken into account when the 
parameter abli\^D is calculated. 


ILLUSTRATIVE EXAMPLES 

Some of the procedures outlined in the preceding section 
for calculating the effective stiffness of supports will be 
illustrated in the solution of two common cover-plate 
stability problems. The first example chosen considers the 
type of restraint offered by the webs of a multiweb structure 
and the second considers the effect of one-sided longitudinal 
stiffeners on plate buckling. 

Buckling of a multiweb structure. — When the webs used 
in a multiweb wing are formed from sheet metal, there is no 
assurance that the deflectional restraint provided to the 
beam covers by the formed channel webs is sufficient to form 
longitudinal nodes along the web lines and thus to force 
buckling of the type denoted as case 6. The subsequent cal- 
culations illustrate a simple procedure that may be used to 
investigate the possible occurrence of buckling in the mode 
denoted as case 3. The calculations apply to a multi web 
beam tested in pure bending and reported in reference 3. 
The beam had four identical channel webs (3 cells) and it is 
assumed that the analysis for a beam with an infinite number 
of cells can be applied. The physical dimensions of the 


beam arc as follows: 

Cover width between webs, h, in 3. 75 

('over thickness, t, in 0. 125 

Channel web depth, bw, in 2. 08 

Channel web thickness, Iw, in 0. 050 

Bend radius between web and attachment flange, in 0. 20 

Diameter of web-attachment rivets, in 3/16 

Pitch of web-attachment rivets, in 9/16 

Distance between midplane of web and line of attachment to 
cover, / (attachment flange assumed to be effectively clamped 
to cover along a line at the inner edge of the rivet shanks when 

closely spaced rivets are used), in 0. 36 

Young’s modulus for the 7075-T6 aluminum alloy, psi . . 10.5X10® 

Poisson’s ratio for the material 0. 333 


111 accordance with the procedure outlined in the preceding 
section, the deflectional stiffness of the channel must be 
evaluated under the action of a sinusoidally distributed 
lateral load of amplitude g on the channel web in the presence 
of the existing bending stresses. This loading is shown in 
figure 14 (a). The lateral loading is applied a distance / from 
the web plane, the distance at which the flange has been 
assumed to be completely fixed to the cover plate. In order 
to compute the deflection at a given cross section, the 
channel is idealized as in figure 14 (b). The attachment 
flange is cut from the web and assumed to be flat and to be 
free of longitudinal compression stress. (This stress is 
usually small in relation to the critical budding stress of that 
portion of the flange between the rivet line and the web.) 
Also, since the buckle length is large compared with the 
distance /, the longitudinal bending stiffness of the flange 
will be neglected in computing the distortions at a given cross 
section. These distortions are shown in figure 14 (c). The 
left-hand edge of the attachment flange is free but maintains 
a zero slope (to match the slope of the attached plate when 
budding occurs in the mode denoted as case 3), whereas the 
right-hand edge is supported against deflection and elastically 
restrained against rotation by the torsional restraint a'. The 
restraint a' represents the resistance to rotation which the 
web offers the flange and is a function of both buckle length 
and the bending stress in the web. Because of the corner 
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radius that~actually exists between the attachment flange 
and the web, the beam cover is assumed to be equally free to 
deflect up or down with the attachment flange. Simple 
tension and compression loading tests on channels with 
corner radii verify this assumption. With these simplifying 
assumptions and boundary conditions, the deflection d at 
any cross section is given by 


&{x) 


l2Dw 


1+4 






1 +' 




^7 


■kx 

T 


The efl'ective stiffness of the channel, defined as the ratio of 
lateral load to deflection, then is 


l2Dy 


l+% 
a / 

1+4% 
« J 


or in nondimensional form 


\j/h^ 

TT+) 



f .-1 

Ow 


J 


e-4 


(7) 


where « is the restraint parameter from figure 9 of reference 
14 and is defined as 


OL ^ hjy 

Djy 


Substitution of the ph 3 ^sical dimensions of the beam into 
equation (7) for 4'b^lir^D gives 


■k^D 0.173e-4 


(8) 


In order to obtain numerical values for the quantit.y 

e must be read from figure 9 of reference 14. Values of e may 
be obtained which are compatible with the bending-stress 
distribution in the beam if the stress in the extreme fiber of 
the web is assumed to be equal to the average stress in the 
beam covers and the lengths of the buckles in the webs and 
covers are equal. From these two conditions, the following 
equations result; 


bw 


b bw- 


= 1.80 


X 

b 


( 10 ) 


The lowest value of the buckling-stress coefficient k which 
simultaneously satisfies equations (8), (9), and (10) is the 
desired value and is found bj' a trial-and-error procedure. 

The first step in this procedure consists in determining by 
trial and error the value of k which satisfies equations (8), 
(9), and (10) for an assumed value of \jh. Values of e are 
read from the curves of figure 9 of reference 14, and values 
of are read from figure 7 of this, report This pro- 

cedure is repeated for several assumed values of X/6. If this 
procedui-e is used, values of k equal to 3.35, 3.25, 3.26, and 
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3.47 are found for assumed values of \/b equal to 0.7, 0.8, 
0.9, and 1.0, respectively. The final step is to minimize k 
with respect to X/6. The minimum value of k for this mode 

of budding (case 3) is thus found to be 3.24 at ^=0.85. 

In order to determine the budding-stress coefficient that 
would be obtained with budding of the type denoted as 
caseV’O" figure 2 of reference 14 may be used to read the 
budding-stress coefficient directly. The use of this direct 
reading chart involves an assumption of an integral joint 
between the webs and the covers, and the indicated k value 
is 4.1, which is considerably higher than the value 3.24 
previously obtained. 

The actual experimental values of the buckling and failure 
stress for the example beam were 

o-cr=33,400 psi 
^failure 36,600 pSl 

and the mode of buckling observed was that of the case 3 
(denoted as wrinkling in ref. 3). If the value Ar=3.24 is 
substituted into the familiar budding equation 

- 


a buckling stress of 34,800 psi is obtained. 

Buckling of a plate with one-sided stiffeners. — In calcu- 
lations of the buckling stress for plates with stiffeners 
attached to one side, the assumption is commonly made that 
the moment of inertia of the stiffeners may be calculated 
about the plane of attachment to the plate. The following 
example illustrates the procedure for obtaining the budding 
stress of the plate-stiffener combination when this assump- 
tion is made and also the slight variation in the procedure 
which is entailed by using the expression from reference 12 
for the effective moment of inertia of a one-sided stiffener. 

Consider the effect of two equally spaced longitudinal 
stiffeners of sturdy cross section on the stability of a long 
compressed plate which is simply supported along the un- 
loaded edges and supported by deflectionally rigid tranverse 
ribs at equal intervals along the length. Assume that the 
stiffeners and ribs offer no torsional restraints to the plate. 


The following physical dimensions are given: 

Plate thickness, t, in. 0. 188 

Plate width between .stiffeners, b, in. 4. 70 

Rib spacing, in. 30 

Cross-sectional area of >4-inch thick Z-stiffener, A„ sq in. 0. 431 

Moment of inertia of stiffener about its centroid, I„ in.'* 0. 203 

Radius of gyration of stiffener, p, in. 0. 686 

Moment of inertia of stiffener about plane of attachment to 

sheet, in.< 0. 524 

Distance between centroid of stiffener and centroid of plate, z, 

in. 0. 956 

Young’s modulus for the 7075-T6 aluminum alloy, psi 10.5X10® 

Poisson’s ratio for the material 0. 333 


The deflectional stiffness of a longitudinal stiffener of 
sturdy cross section is given by (see eq. (2)) 




(»- 


bt by 


if the compressive stress in the plate and stiffener are equal. 
If Eleff is calculated about the plane of attachment of 
El 

stiffener to sheet, is 

12 [1 - (0.33)^] (0.624) E_ 
bD (4.70) (0.1 88)®£’ 

If the buckle length is taken to be the rib spacing, the numer- 
ical expression for \{/PlTr^D is 

/4.70V 0.431 / 30 V,1 

ir*D V 30 / L 4.70X0.188 V4.70/ V 

= 0.1080—0.01195* 

The value of k which satisfies this equation simultaneously 

with the curves of figure 5 at is the desired value. 

® b 4.70 

By trial and error, a common solution is found at * = 3.55. 
In order to verify that * = 3.55 is the lowest buckling-stress 
coefficient, the analysis is repeated by assuming that two 
buckles occur between rib stations. In this particular 
example, this assumption leads to a much higher value of *. 
The buckling-stress coefficient is now computed by assum- 

El 

ing that is given by 


EI,„_Eh 
bD bD 



1 ~\~z. 


4 , 
' bt. 


In order that tlie modal coefficient may be read from 
the curves of figure 13, the buckle length must be assumed. 
The previous calculation indicated that the length of the 
buckle is 30 inches and that it extends across the entire width 

X 30 

of the plate. Thus, with 2 ? = 3, 2=1, and the value 

of Zp, read from figure 13 is 0.80. From the data previously 

EI,,ff . , 

given, IS then 


EI,ff_\2[\ -(0.333)'] (0.203)£' 
bD 4. 70 (0.1 88)3*: 

= (69.4) (2.394) = 166.0 


/0.956V 

V0.686/ 

0.80X0.431 

'*"4.70X0.188. 


With this value for 


ELu 

bD 


the expression 


for \l/b^l-K*D is 


ir^D 


= 0.100-0.01195* 


By the use of figure 5, the value found for * is 3.25. This 
value is about 8 percent lower than the value 3.55 obtained 
when the moment of inertia was rather arbitrarily chosen. 
For other plate-stiffener combinations, the difference in the 
* values calculated by these two procedures can be either 
larger or smaller than the difference obtained in this numer- 
ical example. 
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CONCLUDING REMARKS 

Design charts have been presented which permit the 
evaluation of the compressive buckling stress of a long flat 
rectangular plate with various deflectional and rotational 
elastic line supports running lengthwise of the plate. In 
order to use the charts in a particular plate budding problem, 
the restraint provided by supporting elements such as angle 
and Z-sections and full-depth webs like those used in multi- 
web wing construction must be evaluated. The evaluation 


of the stiffness of these supports has been discussed, and 
possible approaches for obtaining the required stiffnesses are 
presented. Numerical examples have been included to 
illustrate the type of procedures involved in computing 
buckling stresses. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 5, 1953. 


APPENDIX A 


DERIVATION OF STABILITY CRITERIA FOR CASES 1, 2, AND 3 


Although a set of stability criteria could be derived for the 
general case involving any number of lines of support either 
by solving the plate differential equation or by the Rayleigh- 
Ritz energy method, a desirable gain in simplicity is achieved 
by applying the energy method using Lagrangian multipliers 
(see ref. 10) to the individual cases. The latter approach is 
shown in some detail for case 1 , and variations in the method 
are indicated for cases 2 and 3. 

Case 1. — An exact representation of the buckle pattern 
for case 1 is given by the following series 

w=sm-:- 2-1 (Al) 

A 7! = 1,3,5 ^0 

where the origin of the coordinate system lies along a side 
edge of the plate. The sinusoidal deflection along the plate 
center line may be written as 


w(x, h)=A sin ~ 

A 


(A2) 


and the slope along the side edges of the plate may be written 
as 

bw 


/ r> ■ TTX 


(A3) 


Compatibility of equations (Al), (A2), and (A3) requires that 


^ a;, Sin-- — A=0 

71 = 1,3,5 ^ 


y~! a„n—B = 0 

n = l,3,5 


(A4) 


Using equation (Al) permits the so-called strain energj- of 
bending stored in the buckled plate to be written as 

^‘-2 JJ. (X?+V) **= 


P 

4 


Using equation (A2) gives the energy stored in the deflectional 
restraint as 


U2=t [w{x,h)Ydx^t\A^ (A6) 

Jo 4 

and using equation (A3) gives the energy stored in the torsional 


restraints as 


The so-called external work done by the uniform compressive 
load N at buckling is 


‘ 2 Jo Jo \bx) 


dx dy=^- (A8) 

4 ^n = l,‘S,5 


The total potential energ}" may now be written as 

T={U,+U2+Us-V,) 

or 


\h “ 

7’'=4^r= S a. 

TT U 71 = 1,3,5 






(A9) 

where i8=v' and and are the nondimensional 

0 IT U TV U -K U 

deflectional-stiffness and rotational-stiffness parameters, 
respectively. 


The buckling load is determined by the condition that the 
potential energy T' must be a minimum. Since the coeffi- 
cients A and B depend upon the Fourier coefficients, the 
expression to be minimized, is 


Q=T"— A, ( S o,n sin aJ— A s ( S a„n—B^ (AlO) 

\n = l,3,5 ^ / \7? = 1,3,5 / 

where Aj and A 2 arc the Lagrangian multipliers. The 


'' is a minimum when 


dA bB dAi dA 2 

(All) 

y ■ nir 


M ~k — Al sin — — A 2 n =0 


(w=l, 3, 5, . 

. . 00 ) (A12) 


(A13) 

/ 921 A f) 

(A14) 

-= S a»sin^— A=0 

U 7! = 1,3,5 ^ 

(A15) 


«qpL;-- 
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Z1 a„n-B = 0 (A16) 


5 Q 

dA2 n ^173, 5 

Equations (Al2), (A13), and (A14) may be solved for a„, 
A, and B, respectively, and these expressions substituted into 
the compatibility conditions (eqs. (Al5) and (Al6)). This 
substitution results in the following two simultaneous 
homogeneous equations: 


2 A 


sm" 


, mr 


+ s 


A 2 




n sin • 


nir 




= 0 (A17) 


I S 

^ W = ] 


n sin • 


7tTT 




A., 






=0 (A 18 ) 


Each of the infinite sums in equations (A17) and (AI 8 ) arc 
amenable to exact evaluation. Resolving the infinite series 
in equations (,A17) and (AI 8 ) into partial fractions yields the 
following forms: 


sin‘ 


2 2 " 1 


2 “ 1 
7 ^ 

7T 


SDf' 


(I+'t A-‘ «■— * 

\P 4 / 7T 


+ 




S - 


t\1c « = 1 . 3,5 g2 

TT 


n sm 


h = 1,3,5 /I I ^ 

G+t G - 


_2S (-l)-'»sin.^ 


yj. -v/A: 11=1,2,3 


« 2 V 


(- 1 )" 


n sm 


M 

V* —.2.3 ^2 


where 


^ V/3-^/i— 1 


By using equation (6.495) of reference 15, the infinite series 
can be written in closed form. Thus, 


. , nir 

”S,5 


— tan tp—T tanh d 


') 




TV 




tan tanh 6) 


n sm 


mr 


^ 2 _ 7T / 1 

2/3v^ \cos ^ cosh ej 

Substituting the closed forms of the infinite series into equa- 
tions (Al7) and (Al8) and simplifying yield the following 
stability criterion: 


0 = 


2sin <p jsinh 0N 

0 8 — ~ — 

<P I c 


S^k 


tt'^D 


1 

1 ^ 


VCOS <p 

cosh 8/ 


sin (p 

siidi 8\ 

A^|k 

y 



cos <p 

cosh 8/ 

' <A6^ 
t^D 


cos tp cosh d 


2/_J L_^ 

7T \cos (p cosh 8 / 


(A19) 

For given values of k, p, and yblir^D, the value of 
which causes the determinant to vanish is the desired value. 
When the side edges of the plate are simply supported, which 

yh 

is equivalent to setting ~ 2 yT= 0 , the criterion reduces to 

7T IJ 


■k^D' 


Aj_k 

TT^lS 

sin <p sinh 8 
<P , 0 

cos cosh 8 


(A20) 


In reference 16, a stability criterion is presented for the 
compressive buckling of simply supported plates with an 
arbitrary number of longitudinal stiffeners. When equation 
(A7) of reference 16 is applied to an infinitely long plate and 
written in the notation of the present report, it appears as 

4 


1 

jS^V bD 




sm <p 


sinh 8 
8 


cos ——cos (p cos ——cosh 8 
P P 


irq 


which is equivalent to 
4,b^_ 


4y^ 

tt^/3 


x‘*Z) 


sm (p 


sinh 8 

~8~ 


(A21) 


cos ——cos (p cos ——cosh 8 
P V 


I I I 
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when the stiffness of a stiffener is defined by elementary beam 
theory and the stresses in the plate and the stiffener are'equal 
(see eq. (2)). Equation (A21) may be used for plates with 
simply supported side edges and with an arbitrary number of 
longitudinal supports. Equation (A20) may be obtained 
from equation (A21) by substitution of the proper values of 
p and q for case 1 ; that is, p—2 and q=l. 

yb 

For complete fixity of the side edges, 
stability criterion (eq. (Al9)) reduces to 



(A22) 


Solutions of these equations and those to follow are facili- 
tated by a tabulation of the functions tp, sin ip, cos <p, d, sinh 
e, cosh 6 for appropriate values of the parameters k and 
These data are provided in table IV. 

Case 2. — An exact representation of the deflection for case 
2 is given by 

. irx ■ niry , . 

M ;= sm — UmSIU -^- (A23) 

A n = l,2,3 do 


If the same procedure is followed as that for case 1, two 
criteria are obtained, one for symmetrical buckling and one 


TABLE IV 


VALUES OF FUNCTIONS APPEARING IN THE STABILITY 

CRITERIA 


\!b 

<t> 

sin <t> 

cos 0 

e 

sinh 0 

cosh $ 

k = l 

0.4 

6. 08371 

219. 331 

219. 33 

9. 2929 

5430. 37 

5430.37 

.5 

4. 44291 

42. 5051 

42. 517 

7. 6953 

1099. 01 

1099. 01 

.6 

3.31151 

13. 6951 

13.731 

6. 6231 

376. 14 

376. 14 

.7 

2. 45821 

5. 79911 

5. 8846 

5. 8516 

173. 89 

173.90 

.8 

i. 75621 

2. 80881 

2. 9815 

6. 2686 

97. 070 

97. 075 

.9 

1. 10381 

1. 34201 

1. 6736 

4.8115 

61. 454 

61.462 

1.0 

0 

0 

1.0 

4. 4429 

42. 505 

42. 517 

1.2 

1. 1708 

.92106 

. 38941 

3. 8831 

24. 277 

24.298 

1.4 

1.4192 

. 98853 

. 15102 

3. 4764 

16. 156 

16. 187 

1.6 

1. 5209 

. 99875 

. 04987 

3. 1660 

11.835 

11.877 

2.0 

1. 5708 

1.0 

0 

2. 7207 

7. 5626 

7. 6277 

2.5 

1. 5391 

. 99950 

. 03169 

2. 3510 

5. 2004 

5. 2957 

3.0 

1. 4810 

. 99597 

. 08967 

2. 0944 

3. 9987 

4. 1218 

4.0 

1. 3604 

. 97795 

.20885 

1. 7562 

2. 8089 

2. 9815 

6.0 

1. 1708 

. 92106 

. 38941 

1.3853 

1.8729 

2. 1231 

8.0 

1.0390 

.86190 

. .50708 

1. 1781 

1.4702 

1. 7780 

10.0 

.9425 

. 80903 

. 58777 

1.0419 

1.2409 

1.5937 


t = 2 


0.4 

5. 1760i 

88. 4851 

88. 491 

9. 8275 

9268. 36 

9268. 36 

. 5 

3. 40041 

14. 9711 

15. 007 

8. 2094 

1837. 67 

1837. 67 

.6 

2. 0378i 

3. 77171 

3. 9020 

7. 1189 

617.54 

617.54 

.7 

. 449921 

. 465251 

1. 1029 

6. 3310 

280. 86 

280.86 

.8 

1. 4233 

. 98914 

. 14696 

5. 7331 

154. 46 

154. 46 

.9 

1.8232 

. 96831 

-. 24973 

-5. 2619 

96. 422 

96.427 

1.0 

2. 0219 

. 89996 

-.43595 

4. 8813 

65. 897 

65.905 

1.2 

2. 1858 

.81677 

-. 57696 

4. 2995 

36. 825 

36.838 

1.4 

2. 2213 

. 79577 

-.60558 

3. 8737 

24. 050 

24.071 

1.6 

2. 2064 

. 80471 

-.59366 

3. 5467 

17. 335 

17. 364 

2.0 

2. 1240 

. 8.5084 

-.52541 

3. 0735 

10. 786 

10.832 

2.5 

2. 0010 

.90888 

-.41705 

2. 6762 

7. 2305 

7.2993 

3.0 

1. 8857 

.95082 

-.30972 

2. 3977 

5. 4534 

5.5444 

4.0 

1. 6949 

. 99231 

-. 12379 

2. 0264 

3. 7275 

3. 8593 

6.0 

1. 4325 

.99045 

. 13786 

1.6126 

2. 4082 

2.6076 

8.0 

1.2612 

. 95246 

. 30467 

1. 3780 

1. 8574 

2. 1095 

10.0 

1. 1389 

.90817 

. 41859 

1. 2225 

1.5506 

1.8451 


TABLE IV 

VALUES OF FUNCTIONS APPEARING IN THE STABILITY 
CRITERIA — Continued 


X/b 

0 

sin 4> 

cos <t> 

e 

sinh 0 

cosh » 


0.4 

4. 35301 

38. 8501 

38. 862 

10. 219 

13720. 4 

13720. 4 

.5 

2. 2998i 

4. 9359i 

5. 0362 

8. 5831 

2670. 3 

2670. 3 

.6 

1. 0371 

. 86093 

. 50872 

7. 4771 

883. 56 

883. 56 


2. 0685 

. 87868 

-.47741 

6. 6756 

396 42 

396. 42 

.8 

2. 4387 

. 64643 

-. 76297 

6. 0654 

215. .35 

215. 36 

.9 

2. 6095 

. 50733 

-.86174 

5. 5838 

133. 04 

133. 04 

I.O 

2. 6879 

.43829 

-. 89883 

5. 1927 

89. 975 

89. 981 

1.2 

2. 7188 

. 41031 

-.91195 

4. 6934 

49. 410 

49. 421 

1.4 

2. 6786 

. . 44663 

-.89472 

4. 1528 

31. 798 

31.814 

1.6 

2. 6132 

. 50415 

-.86362 

3. 8131 

22. 634 

22. 657 

2.0 

2. 4658 

. 62551 

-. 78021 

3. 3188 

13. 796 

13. 832 

2.5 

2.2932 

. 75022 

66119 

2. 9012 

9. 0705 

9. 1255 

3.0 

2. 1451 

. 83957 

-. 54325 

2. 6067 

6. 7402 

6. 8140 

4.0 

1.9122 

.94229 

-.33481 

2. 2115 

4. 5099 

4. 6195 

6.0 

1.6047 

. 99943 

-. 03390 

1. 7673 

2. 8421 

3. 0129 

8.0 

1. 4081 

. 98679 

. 16198 

1. 5136 

2. 1615 

2. 3816 

10.0 

1. 2692 

. 95486 

. 29666 

1. 3447 

1. 7882 

2. 0488 

j ft=4 

0.4 

3. 51241 

16. 7501 

16. 779 

10. 537 

18864. 3 

18864. 3 

.5 

0 

0 

1.0 

8. 8858 

3614. 34 

3614. 34 

.6 

2. 3416 

. 71735 

-.69671 

7. 77662 

1179. 76 

1179. 76 

.7 

2. 8385 

. 29847 

-. 95442 

6. 9528 

623. 04 

523. 04 

1 .8 

3. 0413 

. 09963 

-. 99502 

6. 3321 

281. 17 

281. 17 

-« 

3. 1221 

. 01949 

-. 99980 

5. 8410 

172. 06 

172. 06 

, 1.0 

3. 1416 

0 

-1.0 

5. 4414 

115. 381 

115. 385 

1.2 

3. 0976 

. 04398 

-.99902 

4. 8273 

62. 433 

62. 441 

1.4 

3. 0106 

. 13062 

-.99143 

4. 3733 

39. 690 

39. 703 

' 1.6 

2 9123 

. 22729 

-.97383 

4. 0240 

27. 953 

27. 971 

' 2.0 

2. 7207 

. 40857 

-.91272 

3. 5124 

16. 750 

16. 799 

2.5 

2. 5133 

. 58776 

-.80903 

3. 0781 

10. 836 

10. 882 

3.0 

2. 3416 

.71735 

69671 

2. 7706 

7. 9528 

8. 0154 

4.0 

2. 0780 

.87411 

48573 

2. 3562 

5. 2280 

5. 3228 

6.0 

1. 7366 

. 98629 

-.16505 

1. 8879 

3. 2270 

3. 3784 

8.0 

1.5209 

. 99876 

. 049876 

1.6191 

2. 4252 

2. 6233 

10.0 

1. 3694 

. 97979 

.20004 

1.4397 

1.9912 

2. 2282 


TABLE IV 

VALUES OF FUNCTIONS APPEARING IN THE STABILITY 
CRITERIA— Concluded 


x/b 

<i> 

sin <i> 

COS ^ 

0 

sinh 0 

cosh 0 

k=b 

0.4 

2. 55191 

6. 37681 

6. 4547 

10. 810 

24768. 18 

24768. 18 

.5 

2. 1587 

. 83210 

-. 55462 

9. 1442 

4680. 05 

4680. 05 

.6 

3. 0604 

.08110 

-.99670 

8. 0123 

1608. 94 

1508. 94 

.7 

3. 3742 

23051 

-.97306 

7. 1881 

6B1. 798 

661. 79S 

.8 

3. 4879 

3.3942 

-.94062 

6. 5580 

352. 43 

352. 43 

.9 

.3. 5123 

-.36227 

-.93206 

6. 0585 

213. 87 

213.87 

1.0 

3. 4928 

-. 34403 

-.93895 

5. 6514 

142. 34 

142. 35 

1.2 

3. 3967 

25235 

-.96763 

5. 0244 

76. 037 

76. 044 

1.4 

3. 2754 

13.341 

-.99105 

4. 5606 

47.815 

47. 826 

1.6 

3.1524 

-.01081 

-.999932 

4. 2010 

33. 369 

33. 384 

2.0 

2. 9270 

. 21295 

-. 97706 

3. 6745 

19. 702 

19. 727 

2.5 

2. 6923 

. 43433 

-. 90075 

3. 2260 

12. 570 

12. 609 

3.0 

2.5019 

. 59695 

— . 80228 

2. 9074 

9. 1272 

9. 1819 

4.0 

2. 2137 

. 80036 

-. 59952 

2. 4767 

5. 9090 

5. 9930 

6.0 

1. 8450 

. 96264 

27078 

1.9881 

3. 5823 

3. 7193 

8.0 

1.6138 

. 99908 

-.042990 

1. 7067 

2. 6646 

2. 8461 

10.0 

1.4520 

. 99259 

. 11852 

1. 5184 

2. 1729 

2. 3920 

/:=6 

0.4 

1. 1164i 

1.3632i 

1.6906 

11.051 

31.542.81 

31542. 81 

.5 

2. 9787 

. 16217 

98676 

0. .3718 

5872. 07 

6872. 07 

.6 

3. 5885 

-.43217 

-.90178 

8. 2285 

1873. 12 

1873. 12 

.7 

3. 7940 

-. 60709 

-. 79461 

7. 3945 

813. 515 

813. 515 

.8 

3. 8468 

-. 64818 

-.76147 

6. 7558 

429. 518 

429. 518 

.9 

3. 8311 

-. 63615 

-. 77155 

6. 2487 

258. 67 

258. 67 

1.0 

3. 7823 

-.59776 

-.80167 

5. 8348 

171.00 

171.00 

1.2 

3. 6459 

-. 48319 

-.87550 

5. 1962 

90. 291 

90.296 

1.4 

3. 4975 

-.34844 

93732 

4. 7227 

66. 232 

56. 241 

1.6 

3. 3548 

-. 21159 

-.97734 

4. 3549 

38. 923 

38. 936 

2.0 

3. 1017 

. 039882 

-. 99920 

3. 8151 

22. 680 

22. 702 

2-5 

2. 8445 

.29274 

-. 95619 

3. 3540 

14. 291 

14. 326 

3.0 

2.6385 

. 48214 

-. 87609 

•3. 0257 

10. 280 

10. 329 

4.0 

2. 3296 

. 72566 

-. 68805 

2. 5808 

6. 50G0 

6.6417 

6.0 

1-9378 

. 93341 

-. 35882 

2. 0745 

3. 9175 

4. 0431 

8.0 

I. 6934 

. 99249 

-. 12230 

1. 7822 

2. 8873 

3. 0556 

10.0 

1.5228 

. 99885 

. 047978 

1. 5863 

2. 3405 

2. 5452 
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for an antisymmetrical wave pattern. Calculations made by considering both modes of budding indicated that, except for a 
very limited combination of values of k and X /6 (k'^A and \jb in the neighborhood of unity), buckling in a symmetrical 
mode requires the highest values of the stiffness parameter to achieve a given buckling-stress coefficient k. Thus, 

for most practical problems, the criterion for symmetrical buckling only need be considered and is given in determinant form; 



sin <p 

1+2 cos <p <p 

I+cos <p 1 

2 “ cos w 


1+2 cosh 9 
1 +cosh 6 


sinh 6 

e 

cosh 6 




3Q^/k 

yb 


/ 

sinh 8 

sin <p \ 

1 ^ 

8 

<p ^ 1 

1 sinh 8 1 
\ 2 

—cosh 6 

sin » 1 1 

2 “ cos ^ 1 


When the plate side edges are simply supported, the criterion reduces to 



(A24) 


Ayfk 

tt"® s in tp sinh 6 

ip 6 


2 — cos <p - — cosh 6 


which is the same as equation (A 21 ) for 5 = 1 , p = 3. 
For complete fixity of the side edges, the criterion is 


(A25) 


4V^ 

xkb^ ^ 


ir^D 


/ 

sinh 8 

sin (p \ 

2 


sinh 8 
8 

1 ^ 

8 

•p V 


sin ip 
‘p 

1 sinh 8 1 

\ 2 

— cosh 8 

sin tp 1 

2 — cos tp 1 

1 

1 

2 ~ cos <p 

i — cosh 8 

^2 1 + 2 cosh 8 

sinh 8 

2 1 + 2 cos tp 

sin tp 

y 



1 +cosh 8 1 
2 

n 

—cosh 8 

^ i + cos tp 1 
2 

— COS ip 


Case 3. For the plate with many lines of support running longitudinally (case 3), the stability will not be influenced 
by the side-edge conditions. Correspondingly, the following function is used to describe the deflection surface: 


W-- 


. TTX 
=sm— 

A -t o j 




(A27) 


wheie the origin of coordinates is taken midway between any two lines of support. Physically the problem thus consid- 
ered is the buckling of an infinitely wide plate column of length X restrained against deflection along continuous longi- 
tudinal lines which are equally spaced across the width of the plate. The stability criterion for this case is 


<kb^ ir^/3 

?r® sin ip sinh 0 (A28) 

<P 8 

1 — cos <p 1 — cosh 8 


which is the same as equation (A 21 ) for <z=l, 00 . 
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APPENDIX B 

DERIVATION OF STABILITY CRITERIA FOR CASES 4, 5, AND 6 


A direct way of obtaining stabilit 3 '^ criteria for cases 4, 5, 
and 6 is by application of the principles of moment distribu- 
tion to the stability of plates as described in reference 7. For 
a long plate supported along longitudinal lines by nondeflect- 
ing supports, the stability criterion is obtained by setting the 
sum of the stiffnesses of the members entering the joint at a 
given support equal to zero. The plate stiffnesses are denoted 
in reference 7 by the s 3 unbol S, with appropriate supemcripts, 
and the carrj^over factors are given by the sj'-mbol C, with 
appi-opriate superscripts. These sj'mbols and their super- 
scripts will be used as defined in reference 7. The support 
torsional stiffnesses a and 7 as defined in this report have an 
absolute value four times as large as S. 

Case 4. — For neutral stability, the sum of the plate stift’- 
nesscs and the support stiffness at the joint along the plate 
centci' line must equal zero. The sum of tbeso stiffnesses is 


~ a+2S' = 0 


(Bl) 


If equation (12) of reference 7 is used, equation (Bl) may be 
wi’itten as 

1 , 25“ 


1-C^ 






which can be put into the following nondimensional form : 

8 

(B2) 


cth 


1 - 


p 




jJl 

t-d 

Solutions to equation (B2) may be readil.v obtained bj^ using 
the tabulated values of S^^blD and T given in reference 11. 
For the particular case of simple support along the plate 
yb 

side edges, equation (B2) reduces to 


ah , 8 




'n + 


D 


= 0 


(B3) 


lb 


With complete fixity of the side edges, “ > equation 

(B2) reduces to 


8 S^^h 


l-(7' 


=0 


or, making use of equation (13) in reference 7, gives 


ah ^ Sh „ 

~D 


(B4) 


With the aid of the tabulated values of and S6/Z) 

given in reference 11, equations (B3) and (B4) have been 

Nb- 

plotted as curves giving the buckling- load coefficient 

7r"/y 


as a function of for constant values of abj-Tr^D. These 
curves are presented as figures 8 and 9. 

Case 5. — If the stiffnesses of the members meeting along 
one of the intermediate lines of support (fig. 2) is summed, 
the following equation for neutral stability is obtained: 


~ a + 5'^ + 5''''=0 


(B5) 


With defined b^' equation (12) of reference 7, equation 
(B5) may be written as 


a-f- 


gu 






1-C' 


^»+y7 


wliich can be written in the nondimensional form 


ab 

■k'^D 


+ - 


4 8^'b 


. _ 4 

O'- ^■k'- D ^ 

'S^^b 


(Bb) 


l-|-4i 

T- lb 

tt’D 


'Phis stabilitx’ criterion is readify solved by using the tabu- 
lated values of S'^blD, S^^b/D, and C given in reference 1 1 . 
yb 

When equal to zero, equation (B6) reduces to 


ab 4: (S^^b s^^b\_ 


yb 

and when 2 jj~ 


», the stabilit,y criterion is 
ah , 4 /S6 , 


(B7) 


(B8) 


Equations (B7) and (B8) have been plotted in figures 10 and 
11 and are presented as curves giving the buckling-load 
Nb^ 

coefficient k= as functions of X/b for constant values of 
TT ry 

ahjic’U. 

Case 6. — For a plate witli many longitudinal lines of 
support (case 6), the condition that the stiffnesses at a joint 
must vanish for neutral stability is given bv 


^ a + 2S''^=0 


(B9) 


In nondimensional form, equation (B9) ma^^ be written as 

(BIO) 


ah , 8 S^''b_^ 


With the aid of the tabulated values of S"'6/Z) given in 
reference 11, equation (BlO) has been plotted as curves giving 

Nb- 

the buckling-load coefficient A-=— ,-y^as a function of X/b foi- 

TT~U 

constant values of ab/ir'-D. These cui'ves ai-e presented in 
figure 12. 
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